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Atmospheric Stirring Measured by Precipitation, 

By Lewis E. Eiohardson. 

(Communicated by Sir Napier Shaw, F.E.S. Eeceived February 26, 1919.) 

Since on the average the water-content of the atmosphere is not increasing, 
the water which descends as precipitation must have been stirred up into 
the atmosphere. So we can compute the eddy-diffusivity from a knowledge 
of the mean precipitation and the mean vertical gradient of mass of water 
per mass of atmosphere. 

Eddy-diffusivity is a statistical measure of the effect of circulatory motions 
which we cannot, or do not wish to, consider in detail. Since the mean 
values in this paper are taken over the whole globe, it follows (unfortunately) 
that even the largest circulatory motions between the pole and the equator 
are not considered in detail, and that their effect is included with the effect 
of the rising currents in cyclone or anti-cyclones, with cumulus-eddies, and 
with the smaller eddies due to dynamical instability, in producing the 
coefficients of eddy-diffusivity at which we shall arrive. 

The General Equations for Stirring. 

In Taylor's theory of atmospheric stirring,* the density of the atmosphere 
and the stirring coefficient are treated as independent of height. 
Under these restrictions he arrives at the equation 

dt ap ' ^^ 

where t is time, h is height, K is the eddy-diffusivity, and x i^^J be either 
the potential temperature, or the vapour pressure, or the horizontal velocity 
in a fixed azimuth. 

The present paper deals with a range of height involving considerable 
variations of density and very large variations in the stirring, so that it is 
necessary to find an equation more general than (1). At the same time, it 
will be convenient to arrange to have a simple expression not only for 
dx/St, but also for the vertical flux. 

Here we may usefully bear in mind the analogy with the conduction of 
heat in a solid. The total water in a portion of air, or the total entropy in 
it, are not altered by gently mixing it ; and the same is true of its horizontal 
momentum in a fixed azimuth, if the mixing takes place without the action 

* 'Phil. Trans.,^ A, vol. 215, p. 3 (1915). 
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of an external horizontal force in this azimutk If, then, Z stands for any one 
of these three quantities or for any other quantity which 

Has its total, for a definite portion of air, unchanged by thel 
internal rearrangement of that portion, J 

Or by delay, (3) 

then the rate of increase of the total Z in the portion must be due to Z 
flowing in over the sides. If we now take the portion to be a horizontal 
layer, and define the " upward flux '' as 

The ratio of amount of Z rising across a large horizontal"!^ 
surface in unit time to the area of this surface, J 

then rate of increase of Z in the layer dh must be 

3 

— ,;--{ upward flux) cl/i. (5) 

If now %, without the dash, be defined as the amount of Z per unit mass of 

atmosphere, then 

The amount of Z in dh is %p(^A', (6) 

so that J^^Ml = — ^ (upward flux). (7) 

vt Oil 

Now we can define a coefficient c such that 

Upward flux = —c^, (8) 

and it will then be for observation to show whether c is a constant, as would 

be nice, or whether c depends on various factors, and in particular on dx/^^- 

But, if e is going to be useful, it must not have to be infinite when 

dx/Sh = 0. Therefore we had better confine ourselves to conditions in 

which 

The upward flux vanishes when 3%/8^ = 0, (9) 

ISTow this implies that 

% must be unchanged by the simple transportation of air"! 
to a different level. J 

Now % is the total Z per unit mass of a definite portion of air. So, as the 
mass is unchanged by transport, (10) implies that 

The total Z in a definite portion of air must be unchanged"! .^-^^ 

by its removal to a different level. J 

Fortunately, (11) is satisfied by the same three quantities which previously 
satisfied (2), so that (10) is satisfied when % is either mass of water per unit 
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of atmosphere^ entropy per unit mass, or horizontal momentum per nnit mass 
(better known as Yelocity), in a fixed azimuth. (11a) 

If, in place of (8), we had made the flux proportional to ^^ d(px)/Bh where 
PX is the total amount of Z per volume instead of per mass of atmosphere 
then, to avoid an infinite value of e% the flux would have to vanish when 
d(px)ldh vanished. That would not be possible for the three given meanings 
of Z, because rising air changes its volume. So it appears that we had 
better keep to the definition (8). 

Inserting (8) in (7), we have 

How the left-hand side of (11b) is equal to p3%/8t+%3p/3^, and by the 
equation of continuity of mass, 

Xdpjdt =: -^XdiiYm, (12) 

where m is the vector momentum per unit volume. 

Let UB now suppose that any time change of mean density of the layer is 
either negligible, as will usually be the case, or else that it is separately 
accounted for. Then (11b) reduces to 

l=;s('l)=-4("5f)' <^'' 

since dp = --gpdh.^ (14) 

Equation (13) is in effect identical with one mentioned by Gr. I. Taylor for 
the case of potential temperature/** 

On the right hand side of (13) there are two variables, p and h, both used 
as measures of height. It would be neater to use only one such measure* 
Therefore, as the dp is indispensable, let us replace the dh from (14). At the 
same time, in order to get a single constant, let us put g^pc:=^^, then 
(13) becooies 

8% __ 3 /^3%\ /-iKx 

3^ dp\ dp/' 

where x ^W have any one of the meanings (11a). Here | is left to be 
determined by observation, but we know beforehand that it is finite when 
3%/3|? = 0. 

The flux is given by (6), (8), and (14) that is to say, if x is the amount of 
Z per unit mass of atmosphere, then the amount of Z rising across a large 
horizontal surface per area per time is 

+ £ ?r^ or, equivalently, — 11- -J^ . (16) 

gdp ^ "^ ' g^pdh ^ ^ 



^ i 
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For the reasons already stated, it appears that it is the quantity ^ of eqtmtions 
(15) and (16) which should he observed, tabulated, and remembered, in preference 
to K of equation (1), which should be regarded only as a stepping stone whereby 
to reach f. The relation between K and f is found by equating the corre- 
sponding expressions for 3%/8^, thus 

diiP'^^A^dpr ^^^ 

If I should happen to be independent of height, equation (17) simplifies to 

t=gyK, . (18) 

and in C.G.S. units the coefficient ^^p^ is numerically equal to unity at an 
average height of 2 kilom. above mean sea level.* 

The physical dimensions of f are 

(Mass)^, (length)"^^, (time)™'l (19) 

Mote on Local Differences of Pressure in a Wind, 

(x. I. Taylorf arrived at equation (1) by considering small portions of air 
which wandered from levels at which they were average samples of the 
atmosphere and arrived at a new level, there to remain by mixing with the 
surroundings. When referring to the diffusion of momentum, he provides a 
proof to show that the local differences of pressure, which Stanton's observa- 
tions have shown to be of considerable magnitude, do not nevertheless on the 
average affect the momentum of a sample during its journey. The local 
differences of pressure, therefore, leave equation (1) unaltered both in form and 
in the value of the eddy viscosity. But Taylor regrets that his proof of this 
applies only to motion in two dimensions. 

In the present paper the difficulty does not arise in the same form. Mo 
attempt is here made to show theoretically that c is numerically the same for 
all three meanings of ^, That question is left to observation. Taylor, who 
originated this stimulating hypothesis, derived it from the consideration of 
" the mean vertical distance through which any portion of the atmosphere is 
raised or lowered while it forms part of an eddy till the time when it breaks off 
from it, and mixes with its surroundings." This distance would be difficult 
to define precisely, as mixing is a gradual process. Consequently, the 
hypothesis seems in need of empirical confirmation, and we come back again 
to observation. 

What has been taken for granted in deriving equation (15), and has now 

* Computer's Handbook, diagram to face Part II, p. 55. 
t 'Phil. Trans./ A, vol. 215 (1915X 
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to be examined more closely^ is that (2), (3), and (9) are satisfied when Z is 
momentum and % velocity, both in the same fixed horizontal direction. Now 
(3) is satisfied because the external forces on a horizontal lamina, due to 
local pressure differences, are normal to its surfaces and therefore cannot 
alter its horizontal momentum. Condition (2) is satisfied because forces 
between the parts cannot alter the momentum of the whole. In examining 
(9) we must distinguish between the mean velocity t?, which here plays the 
part of %, and v\ the deviation from the mean. The question is whether, if 
dv/dh vanish at any level, can there be an upward flux of momentum across 
that level ? In other words can there be an eddy-shearing-stress where there is 
no rate-of -mean-shear ? The converse, of course, occurs where there are no 
eddies. Or to put the question over again in another way: imagine two 
adjacent horizontal layers each say 100 metres thick and having the same 
mean velocity ; would it be possible for the faster moving portions to sort 
themselves out and to flock together into one layer, leaving the slower moving 
portions in the other? If molecules did that sort of thing, the occurrence 
would be one of the exceedingly rare exceptions to the second law of thermo- 
dynamics. We may reasonably expect that there is a corresponding law, on 
an enlarged scale so to speak, applicable to the statistical mechanics of eddy 
motion. 

The only other possibility seems to be that portions of air should wander in 
from quite distant regions, giving a part of the vertical flux of momentum 
proporfcional to d^^/dh^ and a term in either (1) or (15) proportional to d^v/dhK 

Precipitation in the Kilometre next above the Ground, 

To come now from the general to the particular, let Z be mass of water, 
then X is mass of water per mass of atmosphere and, following Hertz, will be 
denoted by fi. The mass of water in the form of vapour rising across one 
horizontal square centimetre per second is given by (16) on the average of a 
large horizontal area. And this must be equal to the average rainfall when 
the horizontal area is extended at that level to cover the whole globe. 
Accordingly we have 

-|- 5j == (mean rate of precipitation). (20) 

g p oh 

The observational data for insertion in this equation, to find |, have been 
extracted from Hann's ' Lehrbuch der Meteorologie,' 3rd edition, to which the 
following paged references refer. According to String, /i in free air at a 
height of 1 kilom. has an average equal to 76 per cent, of its value at sea- 
level (p. 233). The mean value of fju at sea-level has been put equal to 
10*0 X 10"^ which is its saturation value at 14*3° C. and 760 mm. Hg. (p. 221). 
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The temperature 14*3° C. is BoTe's estimate of the mean temperature of the 
earth (p. 143). Consequently the average value of dfjbjdh in the first 
kilometre is — 2-4x 10™^ cm "^ 

Next, the mean precipitation over the earth's surface is about 90 cm. year" -^ 
according to Murray and Fritzche (p. 363). For insertion in the preceding 
equations, this must be expressed as 2*9 x 10""^ grm. cm.~^ sec.~\ This is 
the precipitation at the earth's surface. For comparison with the average 
value of dfj^/dh in the first kilometre, we should use the precipitation at a 
height of half a kilometre. Unfortunately, there is not much instrumental 
evidence concerning change of precipitation with height, as elevated rain- 
gauges are so much disturbed by wind (pp. 319 to 323). However, as rain- 
clouds have an average height (p, 280) well above the height we are 
considering, and as rain, though colder than the air through which it falls, 
can hardly condense any moisbure on to itself, owing to the large latent heat 
of steam (p. 320), and as, further, streaks of rain, descending from clouds, are 
sometimes seen to disappear before reaching the earth (p. 303), the evidence, 
as far as it goes, indicates that the precipitation does not vary much up 
to a height of 500 metres. The mean density at this height is about 
1*2 X lO*'^ grm. cm.~^. Inserting the above mean values in the equation (20), 
we get f = l*4x 10^ cm."^ grm.^ sec.~^ as a mean value for the whole earth 
at a height of half a kilometre. This measure of diffusion is in remarkable 
agreement with the mean deduced by Taylor* from an analysis of tem- 
peratures on the Eiffel Tower ; for Taylor found K = lO'O x 10^ whence 
on the assumption that f did not vary with height, equation (18) would give 
I =: ^2p2 y^ iQ& -». x»45 X 10^ C.G.S. units. It is also of the same order as the 
values of g^p^ (viscosity)/(density) found by Akerblom, Hesselberg, and 
Taylor by analyses of the change of wind with height. Thus Taylor's 
hypothesis, that the same constant applies to both cases, is so far confirmed. 

To revert to the exchanges of moisture, as the precipitation does not vary 
much with height between 500 metres and the ground, it follows from (20) 
that the product of | and dfji/pdh must, on the average, be independent of 
height in this range. This statement will not apply to points within the 
vegetation film, because there transpiration will provide a source of d/m/dt 
which has not been considered in forming the equations. How, just above 
the vegetation film, —d^jdli is sometimes strikingly greater than the mean 
value +2*4x10""^ cm.'"\ which has been taken to hold at a height of 
500 metres. For instance, Hill's results (p. 238) at Allahabad in latitude 
25|^ N. indicate an annual mean value of 9/^/8/^ ^ --66 x 10"^ cm.~^ 
between 1*2 metre and 14*0 metres above the ground. Again, in the metre 

•^ ' Boy. Soc. Proc.,' A, vol. 94 (1918). 
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immediately above a grass mooi% I have observed d^/dh = —1500 x 10""^ cm."^^ 
on an almost calm sunny March afternoon in latitude 49"^ N"., when the 
surface soil was moist. This observation was taken with an aspirated 
wet-and-dry psychrometer. It is, of course, only an instance, perhaps an 
extreme one, but it suggests that the mean value of | in the first metre of 
height may be as low as 1000 or even less, which is interesting, because the 
values of | given by the observation of small floating objects to be published 
presently,* also indicate that f is much less near the earth than at a height 
of 1 kilometre. 

Precipitation at a Height of 8*5 KilometreB, 

Hext, it seems desirable to attempt to estimate the stirring coefficient | 
from the precipitation at great heights; for, although the estimate will 
inevitably be a very uncertain one, it offers the only indication available at 
the moment of the order of this important quantity. The precipitation at 
8-5 kilom. probably does not consist of anything except the slow descent of 
the cloud itself. The evidence for this statement is that the size of the 
precipitating particles is observed to decrease with height, so that, at the 
highest mountain stations, at a height still much below 8*5 kilom., there is 
observed nothing larger than drizzle and fine powdery snow (p. 358).f 
Also that such a cloud-name as " cirro-nimbus '' is not to be found in a list of 
cloud-forms (pp. 269 to 273). We can calculate the velocity Vh with which 
the cloud descends, by Stokes' formula, when we know the size of the ice 
particles. This has been measured by Pernter, by observation of coronae, 
and he has found the diameter to range from O'o x 10"^ cm. to 2'0 x 10"^ cm. 
(p. 257). It will be convenient to make the following calculations separately 
for these two extreme diameters. Now Stokes' formula for the steady 
descent of a very small sphere is 

Y 9 (Diameter)^ {(density of sphere) -—(density of fluid)} 

18 (Molecular viscosity of the fluid) 

The molecular viscosity of the air at the temperature 230° A, which prevails 
at 8*5 kilom., is about 148 x 10""^ C.G-.S. and accordingly, for ice spheres of 
density 0*9, the terminal velocity is Vh= -~0-33x 10"^ x (diameter)^, so 
that, for Pernter's limiting diameters, Vh = —0*08 cm. sec.""^ and --1-3 cm. 
sec."~\ 

In comparison with these velocities due to gravity, those due to the 
electric charge and potential gradient are probably quite negligible, for at a 
height of 8*5 kilom. the mean potential gradient is probably only about 

* In " Weather Prediction by Humerical Process," Cambridge University Press. 
t Paged references here are to Hann's * Me teorologie,' 3rd edition. 
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0*04 volts per centimetre (p. 740), while the mobility of the large atmospheric 
ions has been found by Langevin and Pollock to be less than 10"^ cm.^ 
sec. "^ volt "^ (Kaye and Laby's tables), and the mobility of charged cloud 
particles would probably be even smaller. 

On the other hand, the vertical velocity of the air may sometimes exceed 
the downward velocity of the particles through the air. The vertical 
velocity of the air, or more strictly its product with density, is here treated 
as part of the eddy motion, because its mean value over a level surface, 
including the whole globe, is zero. The air motion, therefore, does not 
affect our measure of precipitation, at least if horizontal variation of 
density can be neglected. 

Suring's table for fju stops at 8 kilom. Taking ^ as 10*0 x 10"^ at the 
earth's surface as before, it follows that mean d^jdh = — 0*3 x 10"^ cm.~"^ 
at 7*6 kilom. (p. 233). However, it is found later that we require dfi/dh at 
the height of the centre of mass of cirrus, cirro-stratus and cirro-cumulus 
clouds, a height which is about 8*5 kilom. Making a short extrapolation 
from Suring's table, it is found that d/jb/dh = —0-2 x 10~^ cm.~Vat 8'5 kilom. 
Even if this figure were ten times in error, the stirring deduced from it would 
still be of interest. 

The density of the air at 8*5 kilom. has a mean value of about 
0*50 x 10"^ grm. cm.~^.* 

The mass of ice per unit, horizontal area of cirrus, cirro-stratus, and 
cumulus clouds has been estimated in another paperf as from 

0*22 X 10"^ to 0*9 x 10""^ gTm.-cmr^ 

according as the particle diameter was 

0*5 X 10"^ or 2*0 x 10""^ cm. 

these being the limits of the range observed by Pernter (Hann, p. 257). It 
will be convenient to follow these limits throughout the argument. The 
corresponding velocities of the particles in still air are, from Stokes' 

formula : — 

Vh= — 0*08 or — 1*3 cm.-sec.""\ 

Now the cloud, when it occurs, often consists of thin layers separated by 
clear air. So, to find a mean mass of ice per volume, we must take, as the 
thickness of the layer, not the thickness of the cloud itself, but the mean 
range within which it occurs, equal to, say, twice the '' probable variation " 
of its height. Now the extreme observed heights for these three kinds of 
cloud are 20 kilom. and 2 kilom. (p. 281). A quantity which varies at 

* Sir W. N. Shaw, * Principia Atmospherica,' Table III. 
t " Measurement of Water in Clouds," infra, p. 19. 
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random according to the law of error, as we may suppose the height of the 
cloud to do within this range, commonly has a " probable variation " eqnal 
to about 1/8 of the total range, which in this case is 18 kilom. So we must 
imagine the ice as spread out in a total range of height of about 4*5 kilom. 
Thus, in the parts of the sky where these clouds are present, the mass of 
ice per volume would be respectively 0*05 x 10"^ and 0'20 x 10"^ grm. cm.""^ 
for Pernter's limiting diameters of particles. Next we require the mean 
fraction of the sky of the whole globe, which is covered by cirrus, cirro- 
stratus, or cirro-cumulus. This fraction will be less than the mean cloudi- 
ness observed at mountain stations at heights of 2 kilom. or 3 kilom., for the 
latter includes also the alto-cumulus and alto-stratus, as well as some 
cumulus and special mountain clouds. 

Some mean annual cloudinesses are set out below :— 



Station. 


Height. 


Latitude. 


Fractional 
cloudiness. 


Beference to Hann. 


Theodul pass ... 
Pic du Midi ..* 
Blue HiU 
Mexican tab] eland 
Peruvian tableland 

- 


km. 
3-3 
2-9 

2-2 
2-4 


46° H. 

42° N. 

42^ N. 

20 -5° N. 

16° S. 


0-64 
0-48 
0-45 
-48 
0-37 


Page 293 
„ 293 
„ 288 
„ 294 
„ 292 



On the average we shall probably not be far wrong if we take 0*3 of the 
globe's sky to be covered by cirrus, cirro-stratus, or cirro-cumulus. It follows 
then that the mean mass of ice per volume of atmosphere at a height of 
8'5 kilom. is respectively 

= 0-015 X 10"^ and 0-060 x 10-« grm. cm.-^ 

for Pernter's limiting diameters of particles. Inserting this, together with 
p = 0*50 X 10""^ grm. cm.~^, 8/i/8A = 0*2 x 10"^ cm.""^ and Vh, respectively 
— 0-08 cm. sec.""^ and —1-3 cm. sec.~\ in the formula which follows from 
equation (20), namely in 



^=fp 



dfijdh 



(21) 



there results | = 3 or 180 cm."^ grm.^ sec.~^ according as we take the smallest 
or the largest ice particles observed by Pernter ; for ^ deduced in this way is. 
proportional to the cube of the diameter of the ice particle. There have 
been many uncertainties in this calculation of f at 8*5 kilom., but even its 
upper limit,* f = 180, is so very much smaller than the value f = 140000 at 
a height of 0*5 kilom. as to make it certain that f diminishes greatly with 

VOL. XCVI. — A. 
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height in the upper* part of the troposphere. In the stratosphere f must of 
course be so small as not to disturb the equality of temperature. 

Sitmmary, 

The equation for diffusion is investigated in the general case in which the 
atmospheric density and the degree of turbulence must both be regarded as 
varying with height, and it is found that the stirring is best measured by the 
coefficient ^ in the equation 

dt dp\ dpi 

in preference to K in 3%/3^ = ^'^xldJi^, 

Here p is pressure, h the height, and % may be either entropy* per unit mass, 
or mass of water per unit mass, or horizontal momentum per unit mass 
(better known as velocity) in a fixed azimuth. 

From estimates of precipitation and of vertical gradient of water per mass 
of atmosphere, as means over the whole globe, it is found that f has the 
following mean values :— 

Height. 
8500 metres 3 to 180 om.""^ grm.^ sec.""^ 

500 ,, 140000 

0*5 „ Possibly as low as 1000, or even less. 

The value at 500 metres is in fair agreement with ^ as deduced from K 
calculated by Akerblom, Hesselberg, and Taylor for heights of a few hundred 
metres. The values of f at other levels are remarkably smaller. 

* Mr. W. H. Dines points out in regard to entropy, that the statement (2) is only 
an approximation, for mixing depends ultimately on conduction or radiation, and these 
are irreversible processes wMch cause the total entropy to increase. However, when 
the initial difference of temperature between the ingredients is a very small fraction 
of the absolute temperature of the mixture, as in the case considered, the error 
committed must be very small. Potential temperature does not involve this error. 



